Reentrance effect in a graphene n-p-n junction coupled to a superconductor 
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We study the interplay of Klein tunneling (= interband tunneling) between n-doped and p- 
doped regions in graphene and Andreev reflection (= electron-hole conversion) at a superconducting 
electrode. The tunneling conductance of an n-p-n junction initially increases upon lowering the 
temperature, while the coherence time of the electron-hole pairs is still less than their lifetime, but 
then drops back again when the coherence time exceeds the lifetime. This reentrance effect, known 
from diffusive conductors and ballistic quantum dots, provides a method to detect phase coherent 
Klein tunneling of electron-hole pairs. 

PACS numbers: 74.45.+C, 73.23.-b, 73.40.Lq, 74.50. +r 



The conductance of a diffusive conductor increases if 
one of the electrodes becomes superconducting upon low- 
ering the temperature, as a consequence of phase co- 
herence between electron and hole excitations (Andreev 
pairs) induced by the proximity to the superconductor. 
The initial increase does not persist to the lowest temper- 
atures. Instead, the normal-state value reappears when 
the thermal coherence length of the Andreev pairs ex- 
ceeds the sample size. This is the so called reentrance 
effect, first observed a decade ago in a metal wire^ and 
in a semiconductor two-dimensional electron gasi^ (The 
theoretical prediction goes back much further££) The 
reentrance effect is now routinely used to measure the 
coherence time of Andreev pairs in a much wider class 
of diffusive conductors, see for example a recent experi- 
ment on multiwall carbon nanotubesi^ It has also been 
predicted to occur in a ballistic quantum dot with prop- 
erly adjusted point contacts^ We refer to Ref. for an 
extensive review of the topic and to Ref. [H] for a tutorial. 

As pointed out by Silvestrov and Efetovjii a p-doped 
region in n-doped graphene (or, converselyly, an n-doped 
region in p-doped graphene) confines carriers in much the 
same way as a quantum dot in a two-dimensional electron 
gas. An electron in the valence band of the p-doped re- 
gion can escape into the conduction band of the adjacent 
n-doped regions, by the process of Klein tunneling. 12 : 13 
This interband tunneling process is highly directional: 
Only electrons near normal incidence are transmitted 
through a smooth n-p interfaced The n-p interface thus 
functions as a constriction in momentum space, analo- 
gously to the constriction in real space formed by a point 
contact in a conventional quantum dot. 

If an n-p junction is in series with a superconducting 
electrode, then the interband tunneling is combined with 
electron-hole conversion (known as Andreev reflection^) 
at the interface with the superconductor. Here we study 
the interplay of these two scattering mechanisms, and 
show that they lead to a reentrance effect in the temper- 
ature dependent conductance. 

We consider a p-doped graphene strip (length L, width 
W), connecting two n-doped regions (see Fig. [1]). One 
n-region is contacted by a normal metal electrode, the 
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FIG. 1: Top panel: Graphene layer connected to one normal 
metal electrode and one superconducting electrode. A bot- 
tom and top gate control the electrostatic doping, so that the 
graphene regions covered by the electrodes are n-type and the 
central region is p-type. Bottom panel: electrostatic potential 
profile across the n-p-n junction. The Fermi level (at Ef = 0) 
lies in the conduction band in the two n-regions, while it lies 
in the valence band in the p-region. 



other by a superconducting electrode. The electrostatic 
potential profile, controlled by a top gate and a bottom 
gate, is assumed to vary smoothly on the scale of the 
Fermi wave length. Close to the Fermi level (chosen at 
Ef = 0) the potential energy increases as U(x) = F(x + 
L/2) near the left n-p interface and decreases as U(x) = 
—F'(x — L/2) near the right p-n interface. The slopes 
F, F' are of order hvkp/d, with v the carrier velocity in 
graphene, fcp the average of the Fermi wave vector at 
the two sides of the interface, and d the thickness of the 
interface. For simplicity, we will take F = F' in the main 
text and consider the effects of two unequal interfaces at 
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the end of the paper. 

The tunnel probabilities T n per mode through a 
smooth p-n junction (kpd 3> 1) have been calculated 
by Cheianov and Fal'koti^ 



T n = exp(-wKvql/F), 



(1) 



in terms of the transverse wave vector q n of the n-th 
mode. These determine the tunnel conductance via the 
Landauer formula, 



Gp-n — 



go}T n , 



go = ±e 2 /h. 



(2) 



The total number of propagating modes in the p-doped 
region is N = kpW/n, with fcp the Fermi wave vector 
in that region. Only a relatively small number N' ~ 
W(kp /d) 1 / 2 of these modes near normal incidence have 
T n close to 1. (In this sense a p-n interface forms a 
constriction in momentum space.) We assume that the 
number N 1 of open scattering channels is still 3> 1, so 
that sums over n may be replaced by integrations over q: 
J2 n ~ * {W/tt) L dq. The resulting tunnel conductance 
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Eq. ([3]) is the conductance of a p-n junction between 
two normal metal contacts. If one of the contacts is su- 
perconducting we can calculate the conductance from the 
formula^ 
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(4) 



Replacing again the sum over modes by an integration, 
this Andreev conductance evaluates to 



Gp_ n — 1.082 G p-n- 



(5) 



The incoherent series conductance Gi nC ohoront of the n-p- 
n junction between a superconducting and normal metal 
contact becomes 
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= 1.040 G . 
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It is slightly larger than the series conductance Go = 
^G p _„ when both contacts are in the normal state. 

To determine the coherent series conductance we as- 
sume that the p-region is weakly disordered (mean free 
path I < L), such that the modes are randomized be- 
fore the electron or hole escapes out of the p-doped re- 
gion. We thus require that the scattering time r = l/v 
is less than the dwell time Tdweii — (L/v)(N/N'), which 
is satisfied if ^Jkpd ^> 1 (assuming kp ~ kp). Dur- 
ing a time r<3 W oii the carrier explores an area LW c f[, 
with W e g = min(W, V^^iwcii) determined by the dif- 
fusion constant D ~ vl. The n-p-n junction corresponds 



statistically to W/W c s quantum dots in parallel, each 
with N e ff = (W e g/W)N' open scattering channels. For 
N c s 3> 1 we may use random-matrix theory to cal- 
culate the average density of transmission eigenvalues 
p{T) through this system. The difference between the 
symplectic ensemble appropriate for Dirac fermions in 
graphene and the orthogonal ensemble of a conventional 
quantum dot does not show up to leading order in N e g, 
so we may ignore it and use the general result in the 
literature for the orthogonal ensemble, 1? 
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(7) 



Combining Eqs. (Q]) and ([7]), the coherent series conduc- 
tance G C oherent evaluates to 



f 2T 2 

Gcoherent = 90 J eff p(T) _ 2 



1.019 G . (8) 



The crossover from coherent to incoherent series addi- 
tion occurs when the thermal energy kT becomes com- 
parable to the Thouless energy 
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with S = 2-k1\v I (LWkp) the mean level spacing (per spin 
and valley) in the p-region. In order of magnitude, Et — 
(hv I ' ^{kpd)^ 1 / 2 ~ h/rdweU- We assume that the gap 
A in the superconducting reservoir is 3> Et- From Eqs. 
© and ([5]) we expect a 2% decrease of the conductance 
upon lowering the temperature T below Ep/k. We will 
now show that this decrease is preceded by an increase, 
such that the conductance is maximal at kT ~ Ep. 

To calculate this reentrance effect we may again use 
random-matrix theory, as in Ref. or we may use 
the equivalent circuit theory of quasiclassical Green's 
functions, as in Refs. HEII Here we follow the latter 
approach. Just as in the random-matrix approach, no 
modification of the conventional circuit theory is needed 
to leading order in N e s. In this quasiclassical limit 
the Green's functions G(r,n) are represented by 4 x 4 
matrices acting in the Keldysh and Nambu spaces, as 
a function of position r and momentum direction n. 
Green's functions Q for Dirac fermions have an additional 
SU(2) x SU(2) structure from the valley and pseudospin 
degree of freedom ; 20 ' 21 which in the quasiclassical limit 
factors out: 



Q(r, n) = hG(r, n) ® (1 + n x ai + n y a 2 ) 



TO- 



(10) 



[The Pauli matrices Tj and a* act on the valley and pseu- 
dospin degree of freedom, respectively, and we take a 
basis in which the Dirac Hamiltonian is v(p ■ <r) ® To.] 

The quasiclassical Green's functions in the supercon- 
ducting lead Gs, in the normal lead Gat, and in the n-p-n 
junction G are 2x2 matrices in the Keldysh spaced 



G = 



R K 
A 



(11) 
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where R, A, and K are 2x2 matrices in the Nambu 
space. In the leads they take their equilibrium values, 
depending only on the excitation energy e: 

R N = ~A N = A 3 , K N = 2[/_ + (1 - /+)A 3 ], (12) 
R S = A S = A 1 , K s = 0. (13) 

Here f± — f(e — eV) ± f(e + eV), with V the applied 
voltage and /(e) = (l + e e / kT )~ 1 the Fermi function. The 
Pauli matrices Ai act in the Nambu space. 

The Green's function in the n-p-n junction has the 
form 

R = Ai cos a + A3 sin a, A = A% cos a* — A3 sin a* , 
K = RF — FA, F = 1 - F+ + F_A 3 . (14) 

The unknown parameters a, F± can be found from the 
conservation law of the matrix current,— 



I N + I s + Is = 

with the definitions 

2g T n [G x ,G} 
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i + T n ({G x ,G}-2Y 



4 = ^[G,A 3 ] 



(16) 

Here [• • • ] is a commutator, {• • • } is an anticommutator, 
and the label X stands for N or S. 

From the retarded component of Eq. (fl"5|) we obtain an 
equation for a, 



tana Z (cos a) = Z(sina) — 
Z(x) = Y J ^T- 1 + x-l)- 1 



2nie 



(17) 
(18) 



The Keldysh component of Eq. ([15)) determines F± , 

/_ Im [cos a Z(sin a)] 



F+ = f+, F- = 



Im [cos a Z(sin a) — sin a Z(cos a) 



(19) 

The electrical current / is related to the Keldysh compo- 
nent (In)k of the matrix current, 



^ / d£Tr[A 3 (7jv)K]- 
2e Jo 



(20) 



The zero temperature differential conductance G(V) = 
dl/dV is then given by 

G(V) = 2 5o (i ;l _//_)coth(Ima)Im[sinaZ(cosa)], 

(21) 

where a is a solution of Eq. (|17p at e = eV and F-/f- 
is given by Eq. (fT9]) . 

The result for G(V) is presented by a solid line in 
Fig. [21 In the limits V — * 00 and V — * we recover, 
respectively, the incoherent limit ([5]) and the coherent 
limit ©■ We find that G(V) takes a maximum value 
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FIG. 2: Differential conductance normalized to its normal- 
state value as a function of the applied voltage, for three val- 
ues of the ratio Gp. n / Gp. n (with Gf. n the tunnel conductance 
of the p-n junction closest to the superconductor and Gp. n 
the tunnel conductance of the other interface). The normal- 
state conductance Go is the series conductance defined by 



( 15 ) l/G = 1/G&, + 1/G| 
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FIG. 3: Same as Fig. [2] but now for the linear response con- 
ductance as a function of temperature. 



G 



1.080 G at eV » E T . Fig. [3 (solid line) shows 
the temperature-dependent linear-response conductance 



G = -2J dsG(e/e)-f(e) 



(22) 



Again we observe nonmonotonic behavior with a maxi- 
mum value G max = 1.058 Go at kT rs Et- 

Figs.[2]and[3]also show results for the more general case 
of two unequal p-n interfaces, for several values of the 
ratio of the tunnel conductances. If the interface closest 
to the superconductor has a higher conductance than the 
other interface, then the reentrance effect is enhanced 
(dashed line), while it is reduced in the opposite case 
(dotted line). This dependence is analogous to that found 
in a conventional quantum dot 4 

In conclusion, we have shown that the conductance of 
a weakly disordered n-p-n junction in graphene coupled 
to a superconductor exhibits a reentrance effect, simi- 
lar to what is found in diffusive conductors and ballis- 
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tic quantum dots: The conductance takes approximately 
the same value at low and high voltages (or tempera- 
tures) and reaches a maximum when eV (or kT) is ap- 
proximately equal to the Thouless energy. This behav- 
ior should be observable in currently available graphcnc- 
superconductor junctions^ and would provide a demon- 
stration of phase-coherent Klein tunneling of electron- 



hole pairs. 
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